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Formation of singularities for a transport 
equation with nonlocal velocity 

By Antonio Cordoba*, Diego Cordoba**, and Marco A. Fontelos*** 

Abstract 

We study a ID transport equation with nonlocal velocity and show the 
formation of singularities in finite time for a generic family of initial data. 
By adding a diffusion term the finite time singularity is prevented and the 
solutions exist globally in time. 

1. Introduction 

In this paper we study the nature of the solutions to the following class 
of equations 

(1.1) 6 t - (HO) X = -vK a 0, xeR 

where HO is the Hilbert transform defined by 

H<)s ip V [m dy , 

ir J x-y 

v is a real positive number, < a < 2 and A a = (—A)~0. 

This equation represents the simplest case of a transport equation with a 
nonlocal velocity and with a viscous term involving powers of the laplacian. It 
is well known that the equivalent equation with a local velocity v = 0, known 
as Burgers equation, may develop shock-type singularities in finite time when 
v = whereas the solutions remain smooth at all times if v > and a = 2. 
Therefore a natural question to pose is whether the solutions to (1.1) become 
singular in finite time or not depending on a and v. In fact this question has 
been previously considered in the literature motivated by the strong analogy 
with some important equations appearing in fluid mechanics, such as the 3D 
Euler incompressible vorticity equation and the Birkhoff-Rott equation mod- 
elling the evolution of a vortex sheet, where a crucial mathematical difficulty 
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lies in the nonlocality of the velocity. Since the fundamental problem concern- 
ing both 3D Euler and Birkhoff-Rott equations is the formation of singularities 
in finite time, the main goal of this paper will be to solve this issue for the 
model (1.1). 

3D Euler equations, in terms of the vorticity vector are 

(1.2) uj t + v ■ Vu = loD(lu) 

where D(u>) is a singular integral operator of uj whose one dimensional analogue 
is the Hilbert transform and the velocity is given by the Biot-Savart formula 
in terms of to. In order to construct lower dimensional models containing some 
of the main features of (1.2), Constantin, Lax and Majda [3] considered the 
scalar equation 

(1.3) u)t + vuj x = ujHlo; 

with v = and showed existence of finite time singularities. The effect of 
adding a viscous dissipation term has been studied in [13], [16], [17], [15] 
and [12]. In order to incorporate the advection term vu x into the model, 
De Gregorio, in [6] and [7], proposed a velocity given by an integral operator 
of uj. If we take an x derivative of (1.1) and define 9 X = uj we obtain a vis- 
cous version of the equation (1.3) with v = —HO which is similar to the one 
proposed in [6] and [7]. 

The analogy of (1.1) with Birkhoff-Rott equations was first established in 
[1] and [10]. These are integrodifferential equations modelling the evolution of 
vortex sheets with surface tension. The system can be written in the form 



z(a, t) — z(a', t) 



< L4) w/^=h pv l 

(1.5) -Tjj; = VKa 

where z(a,t) = x(a,t) + iy(a,t) represents the two dimensional vortex sheet 
parametrized with a, and where k denotes mean curvature. Following [1] we 
substitute, in order to build up the model, the equation (1.4) by its ID analog 

^dT* = ~ H{e) 

where we have identified j(a,t) with 6. In the limit of a = in (1.5) we 
conclude that 7 is constant along trajectories and this fact leads, in the ID 
model, to the equation 

(1.6) e t - (Hd) e x = 0. 

There is now overwhelming evidence that vortex sheets form curvature 
singularities in finite time. This evidence comes back from the classical paper 
by Moore [9] where he studied the Fourier spectrum of z(a, t) and, in partic- 
ular, its asymptotic behavior when the wave number k goes to infinity. His 
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numerical results showed that, up to very high values of k, this asymptotic 
behavior is compatible with the formation of a curvature singularity in finite 
time. Although there has been a very intense activity in order to provide 
a definitive proof of the formation of such a singularity (see discussions and 
references in [9], [2] and [1]) the existing results are mostly supported in nu- 
merics or formal asymptotics and do not constitute a full mathematical proof. 
The same kind of argument was used in [1] in order to show the existence of 
singularities for the ID analog (1.6). 

The system (1.4) and (1.5) with a = has the very interesting property of 
being ill-posed for general initial data. A linear analysis of small perturbations 
of planar sheets leads to catastrophically growing dispersion relations. Several 
attempts at regularization were introduced through the incorporation of effects, 
such as surface tension or viscosity (see [2] for a comprehensive review). In 
the same spirit we will also study the effects of artificial viscosity terms on the 
solutions for our model. More precisely we will prove the existence of blow-up 
in finite time for (1.1) with v = in Section 2 and, conversely, the global 
existence of solutions when v > and 1 < a < 2 in Section 3. 

2. Blow-up for v = 

The local existence of solutions to (1.1) was established in [1]. In this 
section we will show the existence of blowing-up solutions to (1.6) for a generic 
class of initial data. 

Let us consider a symmetric, positive, and C 1+e (R) initial profile 9 = 9q(x) 
such that max x 0o = #o(0) = 1. We will also assume 

Supp(# (z)) C [-L, L] . 

Under these assumptions, it is clear that 9(x,t) will remain positive (given the 
transport character of equation (1.1) for v = 0) and symmetric. Then, H9 will 
be antisymmetric and positive for x > L. This implies the following properties 
for 9(x, t): 

.Supp(0(x,t)) C [— L, L] , 
. max. x 9 = 9(0, t) = 1 , 
■||0|M*) < ||0|MO) , 
.||%»(t)<||*IU>(0). 

Theorem 2.1. Under the conditions stated above for 9q, the solutions of 
(1.1) with v = will always be such that ||0 x ||i,oc blows up in finite time. 

Proof. Since 9 t = -(1 - 9) t = -f t , 9 X = -(1 - 9) x = -f x and H9 = 
—H(l — 9) = —Hf, we can write, from (1.6), 

(2.7) (i-9) t = -H(i-e)(i-e) x . 
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We now divide (2.7) by x 1+5 with < 5 < 1, integrate in [0, L] and obtain 
the following identity: 

r - L (i-6)\ r L {i-e) x H{\ * 



(2.8) 



d_ 

It 



u; 



X 



1+5 



-dx 



7 

J 



-dx . 



Given the fact that 6 vanishes outside the interval [—L,L], we can write 
the right-hand side of (2.8) in the form 



(2.9) 



J o 



(i-e) x H{i-e) 



X 



1+5 



dx 



r 

Jo 



(i-9) x H(i-e) 



„l+<5 



dx . 



In the next lemma we provide an estimate for the right-hand side of (2.9). 



Lemma 2.2. Let f e C^°(1R + ). Then for < 5 < 1 there exists a constant 
C$ such that 



(2.10) 



Jo 



f x (x)(Hf)(x) 



-1+5 



dx > C$ 



1 



o x 



2+5 



f(x)dx 



Proof. First, we recall the following Parseval identity for Mellin trans- 
forms: 



f x (x)(Hf)(x) 1 

,1+5 2tt j . x 



x 1 



A{X)B{\)dX = I , 



with 



A(X)= / 
Jo 

5(A) = / 
J o 



^ A = " 2 fx(x)dx , 



x *-*-i(Hf)(x)dx . 



Integration by parts in A (A) yields 

1 5 s < v 



.4(A) = -(iX- 



2 2' 



x* A 2 *f(x)dx . 



With respect to B{X) we can write 



(A)= / 

7 C 

Jo 

7 C 
7 C 

Jo 



X 2 2 



X lA 2 2 



X 2 2 



X M 2-2 



Ip.v. 

7T 



+ °° /(Q 

■oo X — ^ 

/(0 



de 



f(0 



i r m) i f c 

l /"°° f(f) l Z" 00 f(f) 

71" Jo ^ + ? ^ JO 



d{ 



dx 



_ 7T J X + £ 7T 



x-£ 



dx 



o 



i r 00 x a -2-f 



7T 



X + £ 



dx + 



71" Jo 



00 x iA -i-f 



x-e 



dx 



d£ 

/(0 



dx 



de 



FORMATION OF SINGULARITIES FOR A TRANSPORT EQUATION 



1381 



where we have used Fubini's theorem in order to exchange the order of inte- 
gration in x and £. Using elementary complex variable theory one can write 



1 



, dx + -P.V. 

7T Jo X + i 7T 



00 S**-|-f 



-dx 



■ lim 

_R— »oo 



1 



1 



Z lX " 2 - 



+ 



1 



7T 1 _ e 27«(iA-i-f ) 7^ Z + £ 
1 



7T 1 _ e 2vri(iA^ 



2 2 ^ 



) ^r 3 \{ci,c 2 } 



where T\ and T2 are the paths in the complex plane represented in Fig- 
ures 1 and 2 respectively. Standard pole integration for 1% and the fact that 
Jr 2 \| Cl C2 } = ~~ J/ Cl C2 } i n ^2 (°f- Lemmas 2.2 and 2.3 in [8] where these integrals 
had to be computed for a completely different purpose, for instance) yield then 



h + h 



sin ((-«A + \ + §» 



+ cot 



•x 1 5. 



Hence 



= - 1 + cos((-a + i + f) T ) 

sin((-rt + J + f),r) 




Figure 1: Integration contour T±. 
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and 



27ry_ 0O sm((-a + i + |)7r) 1 2 2 n K > l 



1 

2^2 



oo 

M(A) |F(A)| 2 dA . 



In order to analyze M(A) we define now 

z = a + hi , a = ^- + 7r , 5 = A7r 

which implies, after some straightforward but lengthy computations, 

1 — cosz a sin a + 6 sinh 6 — asinh& + 6 sin a 

2.11 M(A) = z = — + — i. 

sin z cosh o + cos a cosh o + cos a 

Since |-F(A)| 2 is symmetric in A and the imaginary part of M(A) is antisym- 
metric, 



1 f°° 

I= 2^ j /_ oo Re{M(A)}|F(A)| 2 a !A. 



Notice now from (2.11) that 
1 



so that 



(1 + |A|) <Re{M(A)} <C(1 + |A| 



-1 f'OO f'QO 1 

/£ 2^/-J F(A,|2<iA - C U ^ f2(x)dX 



where we have used the Plancherel identity for Mellin transforms: 

,f(x)dx = — / |F(A)| 2 dA. 



/o x 2+&J y ! 2vr 
This completes the proof of the lemma. □ 

Remark 2.3. Inequality (2.10) can be extended by density to the restric- 
tion to M + of any symmetric / G C 1+e (IR) vanishing at the origin. 
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In order to complete our blow-up argument, we have, from Cauchy's in- 
equality, 



o xl+S ~ \Jo x 2+& J \Jo x s 



so that 



1-5 J \J x 2 + 5 



(2.12) j™ ^-dx > C L , 5 ( 

Prom (2.8), (2.10) and (2.12) we deduce 



o x^ 5 



dx 



d f L (1 -9) „ ( f L (1 -6) , 



dt J x l + & ~ ^° V7o x 1+S 
which yields a blow-up for 

Jo x 1+s 

at finite time. Since 

[ L (1-6) r 1-6 [ L dx L 1 - 5 

J -J !^r dx - s T — l x* - T^a s * p| x| 

we conclude that ||^x|Iloo must blow up at finite time. This completes the proof 
of Theorem 2.1. □ 

Remark 2.4. In fact, numerical simulation by Morlet (see [11]) and ad- 
ditional numerical experiments performed by ourselves (see Figures 3 and 4) 
indicate that blow-up occurs at the maximum of 6 and is such that a cusp 
develops at this point in finite time. 

The figures below represent the profiles 6 x (x,t) and 6(x,t) with initial 

data 



6 (x) 

at nine consecutive times. 



(1-x 2 ) 2 , if-l<x<l 
0, otherwise 



3. The effect of viscosity 



Below we study the effect of viscosity {y > 0) on the solutions of (1.1) 
with positive initial datum. First 
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Lemma 3.1. Let 9 be a C 1 solution of (1.1) in < t <T, with a nonneg- 
ative initial datum 9q E H 2 (R). Then, 

(3.13) 1) < 9(x,t) < \\6 \\l°° , 

L 1 7 



2 j. / \\°o\\h 



2u 



(3.14) 2) \\e\\ L i{t)<\\e Q 

(3.15) 3) ||6>||L2(t) < ||6Io||l= and / ||A> 0||f, 2 cft < 

Jo 

Proof. Since 

-| y \A9\ 2 dx = J A9A(H(9)9 x )dx <C\\A9\\ 3 L2 

we have local solvability up to a time T = T(||0o||// 2 (j?)) > (without any 
restriction upon the sign of 9q). Let us also observe that the same result is 
true for the periodic version of (1.1): — ir < x < ir, 



l 2 

We shall prove (3.13) first in the periodic case: Let us define M{t) = 
m&x x 9(x,t), m{t) = min x 9(x,t). It follows from the H. Rademacher theorem 
that the continuous Lipschitz functions M(t), m(t), admit ordinary derivatives 
at almost every point t. Then we may argue as in references [4] and [5] to 
conclude that, at each point of differentiability, M'(t) < and m'[t) > 0, 
implying (3.13). 

Let 4> £ C^(R) be such that <j)>0, 4>{x) = 1 in |ac| < 1 and 4>(x) = when 
\x\ > 2. With R > let us consider 9 Q H (x) = 4>(^)9 Q (x) and let 9 R (x,t) be the 
solution of the periodic problem (1.1) with initial data 9q in — irR < x < ttR, 
< t < T = T(9 ). 

We have that < 9 R (x,t) < \\9q\\l°° with uniform estimates for S/ X 9 R , 
■§l9 R . By compactness, we obtain a sequence 9 Rj , Rj — * oo, converging uni- 
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formly on compact sets to 9, the solution of (1.1) with initial data 6q. Then 
estimate (3.13) follows. 

To obtain inequality (3.14) we proceed as follows: 



^LJedx = J H66 x dx = - j OAOdx = -||A 



ifll|2 

L 2 > 



because f A a 0dx = 0. 
Next, observe that 

\dtj 02dX = J d9 * H6dx - v j 6Aa ° dx 

= ~\j 2 A6dx-v j \A^e\ 2 dx. 

On the other hand 

and the proof of the third part of the lemma follows. 
3.1. Global existence with a > 1. 

Theorem 3.2. Let < 9 e H 2 (R), v > and a > 1. Then there exists 
a constant C, depending only on 6q and v, such that for t > 0: 

(3.16) 1) \\Ah\\ L 2(t)<C , 

(3.17) 2) \\Ae\\ L 2(t) < C(l + t) , 

(3.18) 3) \\A0\\ L 2(t) <Ce ct3 . 

Proof. Integration by parts and the formula for the Hilbert transform 
2H(fH(f)) = (H(f)) 2 -f 2 



yield 

(3.19) 
ld_ 

2dt 



Since 



\kU\ 2 dx = j A66 x H6dx-v J \A2 + i0\ 2 dx 

= - j 9H(e x H9 x )dx - v J \A^ + ^6\ 2 dx 

= ~\j 0(H6 x ) 2 dx + ^ J 6{6 x ) 2 dx-v j \A^9\ 2 dx 

< ||<9o|U°° l|A6»||i 2 — ^ j \A^6\ 2 dx. 



\ko\\ 2 L 2 < R 2 ~ a \\A^e\\ 2 L 2 + R 1 - a \\As + ' e\\ 2 L 2, 
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by taking R sufficiently large and applying inequality (3.15), we obtain the 
desired inequality (3.16). 

Applying A operator to both sides of equation (1.1), multiplying by AO 
and integrating in x, we obtain 

(3.20) ~ J \A9\ 2 dx = J k6k(6 x He)dx-v\\k 1+ ^6\\ 2 L2 

= ~ j\e x ) 2 Mdx-v\\K 1+ ^e\\l 2 

<C j |A6>| 3 dr - z;||A 1+ f #||| 2 , 

where we have used the isometry of Hilbert transform in L 2 , integration by 
parts and finally Cauchy's inequality together with the boundedness of Hilbert 
transform in L 3 . 

In order to estimate ||A#||l3 we make use of Hausdorff- Young's inequality 

(3.21) ||A0|| L 3 < \\M)\\ L% = (^J |^|%)|^)\ 
Picking now a G (1, a) and using Cauchy's inequality we obtain 

2 11 

j_ 1_ 

= T 3 ■ T 3 

For Ii we get 

(3.22) h=f \t\ 2+& \m\ 2 dt+ 1 \t\ 2+ ^\m\ 2 di 

J\£\<R J\€\>R 

<^ 2+s ll^lli 2 + ^l|A 1+ ^ll! 2 . 

With respect to I2 one can estimate 

(3.23) i 2 =! \^\e{i)\di+ [ i^mm 

<[ \0{tM+ f \t\*\om\t\*-*dt 

J\s\<i J \z\>i 

< [ n»dt+ ( [ itf-^dt) 2 ( [ \t\\o(£)\ 2 dt) 2 

J\e\<i \J\e\>i J \J\e\>i I 



'I5l<i \ J \t\>^ ) V\z\>i 

<||e|| L i+Ca||A^||^ 
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From (3.23), (3.14) and (3.16) it follows that 

(3.24) I 2 < C. 
Hence by (3.21), (3.22) and (3.24) we get 

(3.25) ||A0|| L 3 < c^{R^\\e\\\ 2 + -^Wh^ewli) . 

R 3 

To finish let us take R sufficiently large together with (3.20), (3.25) and (3.15) 
to conclude that 

~ J \M?dx<C\%\\l, 

from which (3.17) follows. 
Finally let us consider 

(3.26) IA||A0||2 2= f AdA(6 x He)dx-v\\A 2+ ^d\\ 2 L2 
2 cfct J 

<C||A#|| L =c||A#|| L 2||A#|| L 2 - z/||A 2 +f e\\ 2 L2 

and let us observe that 

(3.27) <c{\\K^e\\l 2 + \\e\\y). 

Therefore, by Holder's inequality, 

||A^||^||A0|| i2 ||A^|| L2 <^||A0||| oo + l||Ae||| 2 ||A^||| 2 , 

and inequality (3.27) we estimate the first term at the right-hand side of (3.26), 
and conclude that choosing 5 small enough, 

|||A0||| 2 <C(||A0||| 2 ||A0||| 2 + ||0||| 2 ) 

which implies the estimate 

||A0||| 2 < ||A0 o ||l 2 e c /o ^h ds + C f \\e\\heK W Ke ^ds. 

Jo 

By (3.15) and (3.17), (3.18) then follows for some large enough C. 

3.2. Small data results for a = 1. In the critical case a = 1 we have the 
following global existence result for small data. 

Theorem 3.3. Let v > 0, a = 1, < 8q £ H 1 and assume that the initial 
data satisfy ||#o||l~ < v. Then there exists a unique solution to (1.1) which 
belongs to H 1 for all time t > 0. 
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Proof. From the previous inequality (3.19) we have for a = 1 

(3.28) ~J \Ah\ 2 dx < ||0 O |U~ \\M\\l> - v J \A6\ 2 dx 

= (||fl ||L--I/)||Afl||i a> 

which implies that if ||0o||l°° < v, then 

r-t 

(3.29) \\A*e\\ L -2{t) < ||A^ ||l= and / \\Ae\\ 2 L2 ds < C\\Ah \\ 2 L 2. 

Jo 

From (3.20) we get 

(3.30) ~ J \A6\ 2 dx< 1 - J \A9\ 3 dx-»\\Ah\\ 2 L2 . 
Since 

||A0||£, < \\Ae\\ 2 L2 ■ \\A6\\ BMO 

and 

||A6>||bmo < C||Af 6»||i, 2 

(we refer to [14] for the corresponding definitions and properties of the functions 
of bounded mean oscillation (BMO)), we obtain 

~ J \A6\ 2 dx<C\\A6\\ 2 L 4A^e\\ L , -v\\A^9\\ 2 L2 

<^\\mh- 

Together with inequalities (3.29) this allows us to complete the proof of the 
theorem. 
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